Survival of the Dirac points in rippled graphene 
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We study the effects of the ripphng of the graphene sheet on the quasiparticle dispersion. This is 
achieved using a generaUzation to the honeycomb lattice of the Momentum Average approximation, 
which is accurate for all coupling strengths and at all energies. We show that even though the 
position of the Dirac points may move and the Fermi speed can be renormalized significantly, 
quasiparticles with very long lifetimes survive near the Dirac points even for very strong couplings. 
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Graphene [l| has been a hot research topic recently [2] , 
primarily due to its Dirac points and the new paradigm of 
rclativistic-like electron dispersion in their vicinity. Such 
dispersion is predicted by almost any hopping model on 
the two-dimensional (2D) honeycomb lattice. However, 
2D systems should not have long-range order [3| , and in- 
deed, free-standing graphene sheets are rippled [J]. In 
epitaxial graphene [5|, coupling to substrate phonons 
becomes important. Both of these can be modeled as 
Holstcin-like coupUng of electrons to out-of-plane opti- 
cal phonons [6|]. Here we investigate the effect of such 
coupling on the electronic dispersion. 

This issue is important because we know from stud- 
ies of Holstein polarons on simple cubic-like lattices that 
even weak-to-moderate electron-phonon coupling has sig- 
nificant effects [7|: while a polaron band with infinite 
Hfctime appears at very low energies, the higher energy 
spectral weight broadens considerably. In other words, 
phonon emission and absorption leads to very short life- 
times for all higher-energy states. This raises the possi- 
bility that the Dirac points, which are at high energies 
above the bottom of the band, may also be "washed out" 
into an incoherent and featureless background. 

We study this problem here, for a single electron. Of 
course, graphene is a half-metal and phonon-mediated 
electron-electron interactions may lead to further broad- 
ening. We assume that such effective interactions, like 
the Coulomb interactions, have little effect on lifetimes. 

Our results show that well-defined Dirac points with 
long-lived quasiparticles are preserved even for extremely 
strong electron-phonon coupling, where most of the rest 
of the spectrum is highly incoherent. Thus, these most 
interesting features are very robust, although their ener- 
gies are shifted somewhat and the slope may be renor- 
malized substantially. These results justify why one can 
ignore the rippling effects and assume a 2D lattice with 
long-lived quasiparticles, as has been done so far. This 
provides a valid description near the Dirac points, but 
would fail if the Fermi energy was anywhere else. 

We use a generalization of the Momentum Average 
(MA) approximation [Sj to calculate the single-electron 
Green's function. MA was shown to be accurate for the 
entire spectrum (not just low-energies) for all coupling 



strengths and in all dimensions, for such problems [9|,llOJ. 
This is so because the MA spectral weight obeys exactly 
a significant number of sum rules, and is accurate for 
all higher order ones. It can also be systematically im- 
proved [10[. We use here the generalizations of the MA^"^ 
and MA*^^^ levels. There is hardly any difference between 
their predictions near the Dirac points, showing that con- 
vergence is reached and we need not go to a higher level. 
The results shown throughout are from MA'^^\ whose 
spectral weight fulfills exactly the first 8 sum rules. 

Consider the honeycomb lattice, with basis vectors 
ai.2 = 3a/2(l, ±l/-\/3), and the three nearest neighbors 
of any site defined by Si = (a,0), ^2,3 = — a/2(l, V3). In 
k-space, the Holstein Hamiltonian for an electron cou- 
pled to an out-of-planc optical phonon mode is: 

n = 5]((/'kc[rfk + h.C.) + n ^(6^6q + St Sq) 
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where c^ and b^{d^ and B^) create electrons, respectively 
phonons on the two sub-lattices. The first term is the ki- 
netic energy of the electron for nearest-neighbor hopping, 
with 0k = —tJ2i=i e'*"^*. Generahzation to other models 
is straightforward. The second term describes the opti- 
cal phonons of frequency fl, for the two sub-lattices. The 
Holstein coupling of the electron to phonons on the same 
site is described by the last term, g being the coupling 
strength. All k-sums are over the Brillouin zone, defined 
by the reciprocal lattice vectors bi_2 = 27r/3a(l,±-\/3). 
We set t = 1, h, ~ I. iV-^oois the number of unit cells. 
Given the bipartite lattice, the single-electron Green's 
function can be defined as a (2 x 2) matrix: 



G(k,uj) 



GH(4 4: 



|0) 



(0|ckG(c.)4|0) (0|ckG(c.)4|0) . 
(0|dkG(u;)c[|0) {0\d^G{uj)dl\0) '' ^ ' 



where G{uj) ~ {cu + iS — H)^^ and |0) is the vacuum. 

The resolvent for the free electron is Go{uj) = {uj + iS~ 
Ho)~^, where Hq = 7i|g=o ■ The free electron propagator 



can be calculated straightforwardly: 

Go(k,a;) = 



Gos(k,c^) e^«WGoA(k,c^) 
e-'«('=)GoA(k,c^) Gos(k,c^) 



where the symmetric and anti-symmetric parts are: 



Gos,A(k,u;) = - 



± 



1 



uj + i6 - \4>{'k)\ w + J5+ |0(k)| 



(2) 



while the phase factor is e'^'^-* = ^k/l^kl- As expected 
in a bipartite lattice, two symmetric bands arise with 
energy dispersions given by ii^k± = i|'/'k| where |(/)k| = 

l + icos'^ (^kyoj + Acqs (^kyo) cos {^k^a) . In 

the first Brillouin zone, the dispersion vanishes at the two 
Dirac points K,K' located at 27r/3a(l,±l/-\/3). 

We now describe briefly the MA'^^) approximation for 
calculating G(k, a;), emphasizing the differences from the 
derivation of Ref. [lO|, due to the two-site basis. Like 
there, first we generate the equations of motion for this 
Green's function, and all the higher order ones it is linked 
to. This is achieved by using repeatedly Dyson's equation 
G(u;) = Go(w) + G(w)T>Go(w), where f = H - Ho is the 
electron-phonon interaction. The first equation is: 



G(k,w) 



i + ;§E^^(k^qi> 



Go(k,c^), (3) 



where Fi (k, qi , uj) is the one-phonon Green's function: 



Fi(k,qi,w): 
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The equation of motion for Fi links it back to G, but 
also to two-phonon Green's functions. The two phonons 
arc both on the same sub-lattice as the electron, or one 
may be on the other sub-lattice. The equations for these 
link them back to -Fi, but also to a new onc-phonon 
Green's function F^ , which has the phonon on the dif- 
ferent sub-lattice than the electron. And, of course, to 
a multitude of three-phonon Green's functions. And so 
on and so forth. All these higher-order Green's functions 
must be proportional to the 2x2 matrix G(k, w) [10|, 
so it is convenient to rescale them accordingly, e.g., 
/i(k, qi,Lj) = g-\/iVG~^(k, Ci;)-Fi(k, qi, w), and similarly 
for all the other ones (see below). Combining this 



J 



/„(qi,...,q„)=.g"A^tG-l(k,c.) 



Ck 



with Eq. ([3]), we find the standard solution G(k, w) 

1 



Go '(k,c^)-I](k, 



where the self-energy is: 



E(k,^) = l^/i(k,qi,o.). 



(4) 



Our task is thus to calculate the momentum average over 
the first Brillouin zone of the generalized Green's func- 
tion /i(k, qi,w). As just discussed, this is the solution 
of an infinite system of coupled equations of motion. We 
cannot solve it exactly, so we proceed to make simplifica- 
tions. At the MA^'^^ level, we replace all free propagators 
that appear in all these equations by their momentum 
averages over the first Brillouin zone. Within MA^^\ we 
keep the equations for /i and f^ unchanged, and make 
the MA approximation only from the second level on, i.e. 
for free propagators of energy co—nD, where n > 2. Either 
of these approximations allows us to solve the resulting 
equations of motion exactly. We proceed to discuss the 
MA'^) solution in more detail. 

As for the simpler case presented at length in Ref. [lO| , 
the simplified MA*^^^ equations of motion can be solved in 
terms of total (and partial) momentum averages of these 
higher Green's functions, over all (all minus one) of their 
phonons' momenta. After such total (partial) averages, 
only contributions from Green's function which have all 
(all except one) of their phonons on the same sub-lattice 
as the electron are non-vanishing. This is because one 
can easily check, using Fourier transforms, that: 



E 

qi,..,qn 



-k-qi-...-q, 
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■<. = ^"^E^''-'''4 
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If any phonon is on the other sub-lattice, the sum over its 
momentum vanishes, since it cannot be at the same site 
as the electron. This shows the variational meaning of 
these approximations [10|jlll|: in MA'"', a phonon cloud 
can appear at any one site. In MA*^^' there can also be an 
additional phonon anywhere else in the system. MA'^' 
thus describes correctly the polaron-|-onc-phonon contin- 
uum [lO[ , but this is a very low-energy feature compared 
to the Dirac points. Diagramatically, both sum all the 
self-energy diagrams, but each diagram is simplified [lO| . 
To summarize, the only Green's functions whose to- 
tal/partial averages remain finite arc (after rescaling): 



G(cL^A.---bl dl_^^Bl...Bl)\0) 



/;^(qi,---,qn)=.9"A^^G-i(k,w)(0| ( ? ] G ( e-'^^-"! dl bLBL.--Bi__ e'^^-^icl „_st.6t ...61 



k-qr qi q2 ■ ■ ■ qn 



k-qT qi q2 ■ ■ ■ qn 



Here, qr ~ X]i=i 1*1 ^-^d for simplicity of notation wc do not write explicitly the k, uj dependence from now on. 



We define the total momentum averages Tn = 
1/A^"^ /n(qi • ■ -qn), and the partial momentum 

averages (5/„(qi) = l/iV""! Eq2...q„ fn{<\\ ■ ■ -fXn) - Tn, 
5/*(qi) = l/A^"-^Eq,...q„/;^(qi---qn). In terms of 
these, the equations of motion are, for any n>1: 

Tn = .90S,n {ng^Tn~l + ^n + l) , (5) 

^/n(qi) =.9os,n [{n- l)g^5/„-i(qi) +(5/„+i(qi)], (6) 
5fn{^i) = 3os,„ [(" - l)9Hr^_Mi) + 5i:,+Mi)] , (7) 
where we use the short-hand notation 

9os,n = gos (^ - nil) = ^ 5Z '^°^ '^'^' ^ ^ '^^) (^) 



where the free propagator is given in Eq. 



cxD, this is a simple 2D integral over the Brillouin zone. 

Recurrence equations of the type ([5])- ([7]) have solutions 
in terms of the continued fractions An{ijj) [1Q |: 



An{uj) 



ngos,i 



1 - g'^90S,nAi+l{^) 



In particular: 



f2=g^A2iL0)fi, 

^/2(qi) = g^Ai{uj^ 

'5/2(qi) = 5'^i(^ 






(9) 



(10) 

(11) 

(12) 



For N 



These can be combined with the exact equations of 
motion for the original (rescalcd) /i(qi) and /*(qi) one- 
phonon Green's functions, which read: 



/i(qi) = G'os(k-qi,cj-f^) 

fnm) = GoAO^-m,oj-n) 



+ ^E/2(qi,q2) +GoA(k-qi,^-r!)l^/*(qi,q2) (13) 

q2 J q2 

+ ^E/2(qi,q2) +Gos(k-qi,o;-r!)^^/^(qi,q2). (14) 



q2 



q2 



r 



From these we can easily derive equations for 
Ti,5fi{qi), Sflim) which combined with Eqs. pO jI - p^ 
allow us to calculate all these quantities. In particular, 
we find Tj ]^,j j^{i) {ll>) = JT^ = 1 • S^mi) (oj), i.e. it is diag- 
onal and momentum independent [12l |. where 



^MAW 



H = 



g'^gosi'^) 



1 - g^gQsi<:i})[A2iLu) - Ai{lu - n)] 



(15) 



with LJ ^ u! — n — g'^Ai{u! — ft). As a result, we obtain 
(5(14,0;) = G'o(k, cj — S^y^^(i)(o;)). This can now easily be 
extended to lattices with even more complex unit cells. 

First, we plot in Fig. [l]the spectral weight A(k, w) — 
— ■^Im{TrG(k, oj)) along high-symmetry cuts in the BZ. 
This corresponds to fairly weak electron-phonon interac- 
tions, with an effective couphng A = g'^/{3tV,) = 0.5. The 
effective coupling is defined such that the crossover from 
large to small polarons is observed for A ~ 1, as usual 
[7[. As expected, the lowest-energy feature is the sharp 
polaron band, whose dispersion flattens out just below 
the polaron-t-one-phonon continuum. Because its quasi- 
particle {qp) weight decreases away from k = 0, it is hard 
to see it in this region. We have verified that all expected 
behavior of the ground-state energy, qp weight, effective 
mass, average number of phonons in the cloud, etc, are 
indeed similar to those expected for Holstein polarons. 

Above the polaron band, we see the polaron-|-one- 
phonon continuum, plus all higher energy features which 



remain centered around the corresponding energies of the 
free electron. Already there is significant broadening at 
all higher energies, except near the Dirac points, which 
continue to be very sharp, indicating long lifetimes. 

This is further analyzed in Fig. [21 where we present 
the imaginary part of the self-energy for two different 
phonon energies, and for different coupling strengths, in 
the vicinity of the Dirac points. Since at the MA^^^ level 
the self-energy is still momentum independent, the char- 
acteristic lifetime is just r ~ — l//m(E). The results in 
Fig. [21 show that as the coupling, and thus the electron- 
phonon scattering, is increased, lifetimes become shorter 
and shorter. However, in a relatively narrow energy in- 
terval, the lifetimes remain very long even at extremely 
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FIG. 1: (color online) Spectral weight along high-symmetry 
cuts in the Brillouin zone, for Q = 0.25i and A = 0.5. 
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FIG. 2: (color online) Imaginary part of the self-energy for 
Q. = Omt (left) and n = 0.25t (right). 



strong couplings which arc well into the small polaron 
regime. In other words, well defined Dirac points still 
exist, although their energy shifts monotonically from 
ujk = at weak coupling to ujn ^ 251 at strong coupling. 
Note that for very strong couplings one can also observe 
the appearance of resonances spaced by Q, which signal 
the Lang-Firsov states, expected when A ^ oo |13l.ll4|. 

Besides the shift in the energy of the Dirac points, we 
also find significant renormalization (decrease) of the ef- 
fective "speed of light" with increased coupling. Also, the 
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FIG. 3: (color online) Spectral weight near the Dirac points 
for a) A = 0.25, b) A = 2.0 and c) A = 4.0 and Q, = 0.07i. 
The red lines indicate the free-electron dispersion (A = 0). 

upper branch has much longer qp lifetimes. These points 



are illustrated in Fig.[3l where we plot the spectral weight 
for k near the K Dirac point at various couplings, to- 
gether with the free electron dispersion (red lines). They 
also agree with the weak-coupling results of Ref. [6( . 

One may wonder if these results are an artifact due to 
the k-independent self-energy produced by MA^^^ jl2| . 
This cannot be the case, since MA^^-* obeys exactly 8 
spectral weight sum-rules for each value of k, so signifi- 
cant spectral weight shifts necessary for the disappear- 
ance of the Dirac points are simply impossible. The 
explanation for this comes from rewriting Eq. ([5]) as 
5os(w) = ^J deg{e)[l/{uj - e + if]) + l/(w + e + iij)], 
where g{e) is the free-electron density of states (DOS). 
The largest contributions come from e ^ ±uj, however 
near the Dirac points ui ^ Q and the DOS vanishes. As 
a result, the self-energy of Eq. (jlSp remains small near 
(just above) the Dirac points, because gos{oj) — > 0. 

We conclude that well-defined Dirac points are a very 
robust feature of the graphene, i.e. rather insensitive to 
rippling effects or, for epitaxial graphene, to the nature 
of the substrate and the buckling due to mismatching. 
This is very fortunate, since it guarantees that the in- 
teresting physics expected because of the Dirac points is 
not affected by such effects. It also explains why they 
can be neglected when studying these quasiparticles, e.g. 
in their interactions with in-plane phonons [15[. 
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